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Abstract 

A primordial degree of circular polarization of the Cosmic Microwave Background is 
not observationally excluded. The hypothesis of primordial dichroism can be quantitatively 
falsified if the plasma is magnetized prior to photon decoupling since the initial V-mode 
polarization affects the evolution of the temperature fluctuations as well as the equations 
for the linear polarization. The observed values of the temperature and polarization angular 
power spectra are used to infer constraints on the amplitude and on the spectral slope of the 
primordial V-mode. Prior to photon decoupling magnetic fields play the role of polarimeters 
insofar as they unveil the circular dichroism by coupling the V-mode power spectrum to the 
remaining brightness perturbations. Conversely, for angular scales ranging between 4 deg 
and 10 deg the joined bounds on the magnitude of circular polarization and on the magnetic 
field intensity suggest that direct limits on the V-mode power spectrum in the range of 0.01 
mK could directly rule out pre- decoupling magnetic fields in the range of 10-100 nG. The 
frequency dependence of the signal is located, for the present purposes, in the GHz range. 
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1 Introduction 



The electron-photon scattering in a magnetized plasma provides computable source terms 
for the evolution of the brightness perturbations [1] which can be studied with diverse initial 
conditions, under various kinds of approximations and in different physical systems ranging 
from the classic problem of line formation of a normal Zeeman triplet [2] to sychrotron 
emission [21 HI El IS] • m conventional Cosmic Microwave Background (CMB) studies the initial 
conditions of the temperature and polarization anisotropies are provided by the standard 
adiabatic mode [71 [5] . The latter requirement implies that the initial radiation field lacks a 
specific degree of linear polarization. By initial radiation field we simply mean, in the present 
context, the initial data for the four Stokes parameters prior to matter-radiation equality 

In the adiabatic scenario, the collision terms of the brightness perturbations for electron- 
photon scattering (see, e.g. [1]) allow for the generation of linear polarization after photon 
decoupling while the circular polarization is, comparatively, not affected, i.e. it is vanish- 
ing both prior to matter-radiation equality and after photon decoupling. The study of the 
linear polarization (and of its potential generation via electron-photon scattering) has a 
long history P [TUl HH [121 H3] dating even before the hypothesis of adiabatic initial condi- 
tions triggered by the formulation of inflationary scenarios. The preliminary detection of 
the polarization autocorrelations [HI [T5] (i.e. the EE power spectrum) as well as of the 
temperature-polarization correlations [TB], [XTJ [TBI E3 Effl EI] (he. the TE power spectrum) 
confirms that the initial data, prior to matter-radiation equality, are predominantly adia- 
batic and lacking any specific degree of linear polarization which arises, to leading order in 
the tight- coupling expansion, because of the quadrupole of the temperature fluctuations. 

In a series of recent papers [221 ES], the idea has been to assume that the initial data of the 
radiation field are dictated by the conventional adiabatic mode but that the electron-photon 
scattering takes place in a magnetized environment. An amount of circular polarization is 
then produced because of the properties of magnetized electron-photon scattering. Analogies 
with the latter phenomenon can be found in the physics of the magnetized sunspots as 
well as in complementary astrophysical situations (see, for instance, 125]). The circular 
polarization arising from unpolarized initial conditions has then been computed and, to 
lowest order in the tight-coupling expansion, is proportional to the monopole of the intensity 
of the radiation field [221 123] (see also [M] for a different kind of derivation). 

The situation explored in [221 [23] does not exhaust the possible sets of initial data for the 
system of brightness perturbations. The statement that the initial radiation field does not 
possess a specified (linear) polarization does not forbid, however, the presence of a primordial 
(circular) polarization. If the pre-decoupling plasma is not magnetized, then, the circular 
polarization will evolve independently both from the temperature fluctuations as well as from 
the linear polarization. Conversely, if the plasma is magnetized, the circular polarization can 
directly affect both the temperature anisotropies as well as the E-mode polarization. 

In the present paper we want to address a situation which is opposite to the one in- 
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vestigated in [2"2"l |2"3] . The idea will be to assume that the initial radiation field has an 
unknown amount of circular polarization. The purpose will then be to constrain the primor- 
dial V-mode by using the magnetized plasma as a polarimeter rather than as a polarizer. 
In [22j 123] the magnetic field acted as a polarizer: the circular polarization was assumed to 
vanish initially and the problem was to compute the resulting V-mode signal. In the present 
paper the magnetic field will act as a polarimeter: the initial circular polarization does not 
vanish and the presence of the magnetic field is used as a diagnostic on the initial radiation 
field. 

The layout of the paper will therefore be the following. In section [2] the relevant governing 
equations will be introduced. In section [3] the V-mode contribution to the temperature 
autocorrelations (TT correlations for short) will be computed and constrained. In section 
H] the same exercise will be repeated in the case of the polarization autocorrelations (EE 
correlations for short). The bounds obtained in sections |3] and @] will be used in section [5] 
with the aim of inferring constraints on the autocorrelations of the circular polarization (VV 
power spectra for short). Section O contains the concluding remarks and the perspectives for 
future studies. 



2 Brightness perturbations in magnetized plasmas 



The results on the evolution of brightness perturbations in magnetized plasmas will be briefly 
recapitulated and the essential notations introduced. The derivation of the equations can be 
found in [22j [23] and will not be repeated here. The evolution equations for the brightness 
perturbations can be written, in Fourier space, a^|: 



AJ + (ik/i + e) Ai = ip' - ik^(p + e' 



Aio + A^b ~ — S P 



-i67eH(l + M 2 )A V l 



|(l-/i 2 K*-^/.(w)(M 2 -l)A 



Ap + (ik/j, + e')A P = 
Ay + (ikfi + e') A v = e'/i|/ c (w)[2A I0 - S P ] - ~iA V i} 



VI 5 



(2.1) 
(2.2) 
(2.3) 



2 The evolution of brightness perturbations appearing in Eqs. (|2.ip . (|2.2[) and (|2.3p is discussed in the 
conformally Newtonian gauge where the perturbed entries of the geometry read <5 s <?oo = 2a 2 </> and S s gij — 
2a 2 ip5ij. It should be noticed that the conventions are different from Ref. pJ|B] where <f> and tp are exchanged 
and where the signature of the metric is mostly minus. 
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where the prime denotes a derivation with respect to r, i.e. the conformal time coordinat^f). 
Moreover, Sp, f e (u>) and e' are defined as 

S P = Ap 2 + A P0 + A I2 , f e (co) = ^, (2.4) 

UJ 

.a / \ r° ~ a i 

e=x e n e a^ e —, e(r, r ) = / x e n e a 7e — dr. (2.5) 

OjQ Jt CLq 

Concerning Eqs. ( 12.4D and (12. 5 ft few comments are in order. The limit / e (u;) — >■ corresponds 
to the standard situation where the plasma is not magnetized: indeed f e (uj) denotes the ratio 
between the Larmor frequency of the electrons and the angular frequency of the observational 
channel. As usual e' is the differential optical depth while e(r, r ) is the optical depth. It 
is appropriate to remind the explicit form of f c (uj) which does depend upon the angular 
frequency u = 2itv and upon the magnetic field intensity B u 



„. , ur p 10 / B u \ /GHz\ . , e\B ■ n\ 
f (u) = -52 = 2.79 x 1(T 12 + 1 , w Be = J 2.6 

where z* is the redshift to last scattering, i.e. z* = 1090.79lo;g2 according to the WMAP- 
7yr data. In the present paper the values of the cosmological parameters will be taken to 
coincide, without loss of generality, with the ACDM best fit to the WMAP 7-yr data alone; 
the critical fractions of baryons, CDM particles and dark energy will be, given respectively, 
by 

(fi b , n c , fi dc ) = (0.0449 ± 0.0028, 0.222 ± 0.026, 0.734 ± 0.029), (2.7) 

while the rescaled Hubble parameter, the adiabatic spectral index and the optical depth to 
reionization are 

{ho, n s , e re ) = (0.710 ± 0.025, 0.963 ± 0.014, 0.088 ± 0.015). (2.8) 

The power spectrum of the curvature perturbations will be assigned, as usual, at the pivot 
scale k p = 0.002 Mpc" 1 as 

V n {k) = Ak^yY ' ^= (2.43 ±0.11) x 10- 9 . (2.9) 

The system of equations given in (12. ip . f)2.2p and (12. 3 j) can be studied in different physical 
limits. If we assume that the initial V-mode polarization vanishes exactly we can keep, at 
least initially Ayi = and Eqs. (I2.ip - fl2.3p reduce to 



Aio + fWb — S P 



(2.10) 



Aj + (ikfi + e')Aj = ijj' - iky4 + e' 

Ap±^/i± e , )A P = l(l-ti 2 )dSr, (2.11) 
A' Y + (ikp + e')A v = e , /i/ B (w)[2A M) -S , p]. (2.12) 

3 In this paper conformally flat geometries will be considered where the background metric = a 2 (r) 
is written in terms of the Minkowski metric rj^ = diag(+l, — f , — f , — f ) and of the scale factor a(r). This 
choice is directly dictated by the adopted set of fiducial parameters (see below Eqs. (|2.7j) and (|2.8[) ). 
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The system of Eqs. fl2.10l) - fl2.12p describes the situation where the radiation field does not 
have (initially) any circular polarization. As a result of the electron-photon scattering in a 
magnetized environnment a tiny amount of circular polarization is produced according to 
Eq. (I2.12p and has been computed in [221 123] . In the present paper the logic will be opposite 
to the one leading to Eqs. fl2.10p - fl2.12p . More specifically, it will be assumed that the initial 
radiation field is circularly polarized (i.e. Ayi 7^ initially). The question will then be how 
large could Ayi not to affect the temperature and polarization autocorrelations. In [23] the 
coupling to Avi appearing in the analog of Eqs. f)2.2p and (12. 3 p has been incorrectly squared. 
We take here the chance of correcting this coupling which does not affect the results of [23] 
since, in that context, A V i did vanish prior to matter radiation equality. 

The maximum of the microwave background arises today for typical photon energies of 
the order of 10~ 3 eV corresponding to a typical wavelength of the mm. At the time of 
photon decoupling (i.e. z ~ z* ~ 1090) the wavelength of the radiation was of the order 
of 10 _3 mm ~ fim. Since the magnetic field we are interested in is inhomogeneous on a 
scale comparable with the Hubble radius (at the corresponding epoch) the guiding centre 
approximation [27] can be safely employed as already discussed in [22], [23] . There are different 
ways of introducing the guiding centre approximation and the simplest one is to think of a 
gradient expansion of the background magnetic field, i.e. denoting with B the (comoving) 
magnetic field intensity we can write that 

Bi{x, r) ~ Bi(x , t) + (x j - 4)djBi + ... (2.13) 

where the ellipses stand for the higher orders in the gradients leading, both, to curvature 
and drift corrections which will be neglected in this investigation as they were neglected in 
[2"2"| |2"3"]. Higher order in the gradients demand the inclusion of higher multipoles of the field 
(see e. g. [28]). The latter effects can be neglected to lowest order in the guiding centre 
approximation. 



3 Temperature autocorrelations 

The temperature fluctuations can be written as 

A T (n, t ) = °S1 Yt m (n), n = (tf, (p) (3.1) 



£m 

From the line of sight solution of Eq. ( 12. ip . the power spectrum of the temperature correla- 
tions receives two separated contributions stemming, respectively, from the intensity of the 

T (vh 

a \m = a £m + a £m- (3-2) 



radiation field (denoted as a^J and from the circular polarization (denoted as 



To make the notations clear we want to stress that a^2 denotes the V-mode contribution to 
the temperature correlation and not the power spectrum of the V-mode itself (which will be 
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introduced later on in section E]). According to Eq. ( 12.1 p the two terms at the right hand 
side can be written as 

= J^ym I Sk £ *P jf d<pYl m (fi,ip) jH e-^e-^)M(fc,/i,r)rfr, (3.3) 



77 {v) 



(2tt) 3 / 2 



where \i = cos??; the two generalized sources Afi(k,fi,r) and J\fy(k, f/,, r) are given, respec- 
tively, by: 



J\fi(k, fi, t) = if)' — ikfi(p + e' 



Aio + ^v h - ^P 2 (fi)S P 



(3.5) 



Nv(k, fi, r) = ~i e'f e (u)(l + /i 2 ) A V1 . (3.6) 

For angular scales larger than 1 deg, the visibility function can be approximated with a 
sufficiently thin Q Gaussian profile [29], [30J [31] : 

'2 1 



/C( T ) = e'e- e(T ' T0 > = \ e 2 < , (3.7) 

7T (T * 

i.e. To ^> r* and to ^> a*. In the latter limit and assuming that the V-mode and the adiabatic 
contribution are uncorrelated the angular power spectrum at large scales is given by 

Q (TT Wf+^ V) , (3.8) 

where the two terms at the right hand side are, from Eq. (13. 2p . 



Jf ~E <I«2J 2 >, ^ (v) = ^tI (l«Sl 2 )- (3-9) 

m=—l rn=—£ 

According to Eq. (13. SJ) , the temperature autocorrelations contain two physically different 
terms: the first contribution is given by the fluctuations of the intensity of the radiation field; 
the second contribution stems from the V-mode. The second contribution vanishes either 
when the magnetic field is zero or in the case when the circular polarization of primordial 
origin vanishes. To derive Eqs. ( 13. 8 p and ( 13. 9 p it has been assumed that the V-mode 
contribution is not correlated with the adiabatic mode which is instead responsible for the 
perturbation in the intensity of the radiation field |S] . 



4 Note that tq denotes the angular diameter distance to t*. It is appropriate to mention that the visibility 
function has also a second (smaller) peak which arises because the Universe is reionized at late times. The 
reionization peak affects the overall amplitude of the CMB anisotropics and polarization. It also affects 
the peak structure of the linear polarization. The fiducial set of parameters adopted in Eqs. (|2.7I) - (|2.8I) 
suggest that the typical redshift for reionization is z r eion = 10.5 ± 1.2 and the corresponding optical depth 
is e rc = 0.088 ± 0.015. For semi-analytic purposes the second peak can be modeled with a second Gaussian 
profile [32] . 
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It is now useful to derive the evolution equations for the lowest multipoles of the V-mode 
polarization. Recalling the present conventions for the Rayleigh expansion 

A,(k, (a,t) = j £ HY (2* + 1) W \e(k, r), (3.10) 

i 

Ay(k, /x, r) = £ H)' {2i + l)P t (ji) A ve (k, r), (3.11) 

e 

the following set of relations can be easily derived for the monopole, for the dipole and for 
the higher multipoles of the brightness perturbation associated with the V-mode: 

A^ + e'Ayo + kA vl = 0, (3.12) 
A' Y1 + hA Y2 - ^A vo = -jje'Ayi, (3.13) 

A^ + e'A w = ^[Mv^i, -(£+ l)A v(m) ]. (3.14) 

In analogy with what customarily done in the case of the brightness perturbations of the 
intensity it is practical to define 

6y — - [ d/j, d<p Ay (k,u, t) = 4 A V o, (3.15) 

71 J 

8y = — / d/j dip/j Ay(k, u, t) = 3A;A V i, (3.16) 

47T J 

ay = -— J dndip(ji 2 - -^jAy(k,n,r) = 2A V2 - (3.17) 

In terms of Sy, By and ay Eqs. A3. 1 2|) and (13 . 13[) can be written, over large scales, as 

5v + ^v = 0, (3.18) 

9 y + k 2 ay = ^Sy. (3.19) 

It is therefore possible to set initial conditions by requiring that, prior to matter-radiation 
equality, the power spectrum of the dipole does not vanish and it is given as 

Vy(k)=Ayi^—j , k p = 0.002 Mpc -1 , (3.20) 

where k p is the same pivot scale used to assign the adiabatic mode in Eq. (I2.9p . By combining 
Eqs. (I3.12p and ( 13 . 1 3j) the evolution of the dipole can be written as: 

0y + ^-dy = ~k 2 ay. (3.21) 

Setting to zero all the multipoles I > 2, the solution of Eq. ( 13 .2 1 j) prior to equality will be 
given by: 

9y(k, t) = Ci{k) cos (c sv kr) + C 2 (k) sin (c sv kr), (3.22) 
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where c sv ~ l/y/3. Equation ( I3.22p stipulates that there are two separate classes of initial 
conditions: initial conditions where, prior to equality, the degree of circular polarization 
vanishes (i.e. C\{k) = implying that 6(k, r) — > for r — > 0); initial conditions where, prior 
to equality, the degree of circular polarization goes to a constant (i.e. C^/c) = implying 
that 6{k, t) d(k) for r ->• 0). 

The purpose of the present investigation is to set a bound on the amount of circular 
polarization present prior to equality and it is therefore appropriate to postulate that C*2(^) = 
while the constant C\(k) is directly related to the power spectrum of Eq. ( I3.20p . More 
specifically we shall have, for r <C r eq , that |Avi(fc,r)| 2 = 2n 2 Vv(k)/k 3 ; however, because 
of Eq. f)3.16p . and in the same limit, the following chain of equalities holds 

\6 Y (k,r)\ 2 = \d(k)\ 2 = 9k 2 \A Y1 (k)\ 2 = ^Vy(k). (3.23) 

For large angular scales, the sudden decoupling limit implies that the angular power spectrum 
of the circular polarization is given by: 

^ (V) = ^/>)(^y W MM*,*), (3.24) 



2 Vfcp 



' c 

X v (£,n v ) = A 2 (£)X 1 (£,n v ) + B 2 (£)X 2 (£,n v ) + C 2 (£)X 3 (£,n v ) 
- 2A(£)B(£)X 4 (£, n v ) - 2A(£)C(£)X 5 (£, n v ) 

+ 2C(£)B(£)X 6 (£ 1 n v ). (3.25) 

The coefficients A(£), B(£) and C(£) appearing in Eq. (I3.25P are the same as the ones 

obtained in Eq. ( 1A.8I) of appendix lAl The integrals Xi(£,n v ) (with i — 1, ...6) can be 
computed analytically and the result is 

1 r (f-TW £ +lf- 2") 

X 1 (£,n v ) = —=-) f-) f, (3.26) 

l r (f-if) r (^ + if-f 

l 2 (£,n v ) = f-) (3.27) 



l r f-f r £ + f + | 

X 3 (£,n„) = —=-) f-) f, (3.28) 

i K-2)r(|-f)r(£-| + f 

Zi&n.) = \ / x-^, (3.29) 

4v^ r(3 - ^)r(| + 



k-2)r|5-|rf + | + f 



X 5 (£,n v ) = —7-7= 7 — x / x -, (3.30) 

4v?r r3-^m+f-^ 



8 



2 r(§-2*)r(2* + *-i 

X 6 (£,n„) = -= ^ '-A / (3.31; 

k - 6)K - 4)r(j + i - 

where we assumed — 3 < n v < 3. By demanding that the intensity power spectrum is always 
smaller than the corresponding V-mode contribution, Eqs. (13. 8p and (13.91) imply that 

Jf } <jf ^Cf T) . (3.32) 

Equation (13.321) holds in the standard ACDM scenario and under the assumption that the 
V-mode contribution and the adiabatic contribution are not correlated. The inclusion of 
the cross correlation may entail the addition of a further power spectrum and will not be 
explicitly discussed in this paper. From Eq. (13. 3p . after integrating once by parts the term 
—ik/Mj), it is easy to obtain the following expression for af' m 

_(I) _ (~')^m0 

a£m ~ (2tt) 3 / 2 



VW+iVAtt J d 3 kA u (k,T ). (3.33) 



The term Au(k, r ) can be read off directly from Eq. (13.101) and by bearing in mind the line 
of sight solution of the heat transfer equation. The result is 



A M (fc,r ) 



Ai (k,n) + (f>(k,n) 



3i[k(r - n)} + / (0' + V') jt[k(T Q - t)} dr, (3.34) 



where je(x) are the spherical Bessel functions [331 EI]; in Eq. (I3.34[) the first term corresponds 
to the ordinary Sachs- Wolfe contribution while the second term corresponds to the integrated 
Sachs- Wolfe effect which operates, in the ACDM scenario, after matter-radiation equality. 
Over sufficiently large angular scales the difference between (f) and ip can be ignored and Eq. 
(I3.33P can be written as: 



A K (fc,r ) = -^P- ji[k(r - n)} -2 jT (^) K.(k) ji[k(r - r)]dr, (3.35) 

where 



T(t) = \--^ j a 2 {r')dr\ (3.36) 



a 2 Jo 

as it can be easily deduced, by direct integration with respect to the conformal time coordi- 
nate, from the definition of curvature perturbations in the conformally Newtonian gauge: 

K = -tp- n2 _ n , . (3-37) 

and under the assumption, already mentioned, that = ip. It is appropriate to recall that 
the first thorough analytical discussion of the integrated Sachs- Wolfe effect in the ACDM 
paradigm dates back to the results of Ref. [35]. The integral over the conformal time 
appearing in Eq. ( I3.35|) can be performed by recalling the the first zero of the spherical 
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Bessel function is given by k(r — r k ) ~ (2£ + l)/2, i.e. r k = r — (2£ + 1) /(2k). Thus we can 
also write that 

r°fdT\ l/dT\ r 00 

I {^) n *^^ T °- T ^ dT -k[^) T _ U *W J, ^ dx - ( 3 - 38 ) 
Equation (I3.35|) becomes then: 



A u (k,r ) = _^*)^[fc( To - T .)] - (df) TZ^k)^-^^. (3.39) 

5 V (IT / r=Tfc ft p / £+2 



The exact solution of the Einstein equations in a ACDM cosmology can be written, in the 
cosmic time coordinate, as 



a(t) = aijsinh 



-\JVtM)H (t — 1 



2 ' 3 , (3.40) 



The conformal time coordinate r(a) and T(a) (the transfer function) can then be expressed 
in terms of integrals over the (normalized) scale factor a as: 

(\ 2 ^ TV x i V a 3 + 1 ,q ^3/2 

r(a) = ^o7n^i TFtT' T(a) = 1 -^7^io TFtt^- (3 - 41) 

The two integrals of Eq. (13.411) can be performed analytically in terms of elliptic func- 
tions; however the obtained result must anyway be integrated (numerically) over k. This 
will mean that, at the end we will not have a closed analytic expression. We will proba- 
bly gain in accuracy but the simplicity of the results will be partially lost. It is useful to 
notice, in this respect, that the ordinary and the integrated Sachs- Wolfe contributions are 
reasonably well separated in scales. The Sachs- Wolfe contribution typically peaks for co- 
moving wavenumbers k ~ 0.0002 Mpc -1 while the integrated Sachs- Wolfe effect contributes 
between k min = 0.001 Mpc" 1 and A; max = 0.01 Mpc" 1 . For comparison recall that the pivot 
scale k p = 0.002 Mpc" 1 corresponds, for the best fit parameters of the WMAP 7yr alone 
(and in the light of the ACDM scenario to £ p ~ 30. The latter remark suggest that if the 
integrated Sachs- Wolfe contribution is neglected, the ordinary Sachs- Wolfe term leads to a 
more constraining bound simply because it does not always dominate at large scales. The 
strategy will then be to derive an analytic form of the bound by considering only the ordinary 
Sachs- Wolfe contribution. More refined treatments could certainly improve on this aspect 
but, for the present purposes, the estimate will be sufficiently accurate as we shall see in a 
moment. By inserting Eq. ( I3.39|) into Eq. H3.33|) the angular power spectrum becomes 

rd) _ ^(koY^ 1 , „_ 2er 



F? 3 = -^[f) AnG X ^ n s), (3.42) 

i r(§-f)r(Wf-i 

l n (£,n s ) = T -=-) (—? f> ( 3 - 43 ) 



2v^F 



T[2-'f)T[£ + l-f 
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where e re accounts for the overall suppression inherited from the reionization peak of the 
visibility function [32]. In the range 2 < t < 40 the ratio Tn{£, n s )/Xy(£, n v ) can be usefully 




Figure 1: The results of Eqs. (I3.44p and ( I3.47[) are illustrated for different values of the 
spectral indices. 



approximated in a factorized form as 

t id \ - — : i i I — 

r(E,n v ,n s ) - 



2v(*,n.) ( n 2 _ i 2riv + 39)1^2 - f Wf - ^) 



. (3.44) 



The accuracy of Eq. (I3.44p in the physical range of spectral indices is of the order of few 
percent depending upon the value of the spectral indices and upon the specific multipole. The 
bound of Eq. (I3.32p can now be enforced. By requiring, within our fiducial set of parameters, 
that the V^-mode contribution does not exceed the observed temperature autocorrelation we 
are led to the following condition 

where the term A/tt is given by: 

A/"tt = 1.156 x 10 6 x (0.0354)"" e" 2 ^ r(£, n v , n s ). (3.46) 

In Eq. (13.451) Da denotes the (comoving) angular diameter distance to last scattering while 
denotes the redshift to the last scattering. All the typical values appearing in Eqs. (13.451) 
and (I3.46P refer to the WMAP 7yr data alone. It is often convenient, for sake of simplicity, to 
average r(£, n v ,n s ) over a suitable range of multipoles and the resulting averaged expression 
will be 

(/n v -n s +l _ pn v -n s +l-\ 4 T ( | — ^rf J T [ 4 — ^ 



r(n v , n s ) = 1 max " mm ^ '-L^ — ^-j -. (3.47) 

(W - Cin)K -n s + l) ( n 2 _ i 2rlv + 39)r( 2 -| r 



11 



The approximations obtained in Eqs. fl3.44j) and (I3.47P are illustrated in Fig. [TJ where the 
exact results (triangles) are compared with the approximate result of Eq. H3.44[) (open stars) 
and with the average values of Eq. (I3.47[) (open boxes) deduced from Eq. (I3.47j) . The scalar 
spectral index n s has been taken to coincide with the best fit parameter of Eq. (12. 8p while 
the spectral index n v has been illustrated in two specific cases, i.e. n v > n s (plot at the 
left) and n v < n s (plot at the right). The dependence upon the multipole £ can therefore 
be safely accounted for by using the derived approximate expressions and this trick will 
permit further simplifications in the final expressions. The current bounds on the V-mode 
polarization are rather loose for the standards we are used to in the case of some pivotal 
cosmological parameter such as the ones appearing in Eqs. (12. 7p . ( 12. 8 p and (12. 9p . This 
discussion will be approached in section 



4 E-mode autocorrelations 

In the heat transfer equation for Ap, there is a dependence upon Ayi (see the right hand side 
of Eq. (12. 2p ). This observation can be used to infer a further constraint on the primordial 
power spectrum of the circular polarization. For this purpose the brightness perturbation 
Ap(/c, fM, t) must be appropriately related to the E-mode polarization. In full analogy with 
what has been done with the TT correlations in section [3] the idea is to compute sepa- 
rately the standard contribution and the V-mode contribution. The orthogonal combina- 
tions A±(n,r) = Aq(n,r) ± iAu(n,r) transform as functions of spin-weight ±2 and can 
therefore be expanded in terms of spin ±2 spherical harmonics [36j EZ] (see also [38j [391 SO] 
for a background on spin- weighted spherical harmonics): 

A±(n,r) = ^2a ± 2,e m ±2Y tm {n). (4.1) 

im 

The E- and B-modes are, up to a sign, the real and the imaginary parts of a±2/ m , i.e. 

a £m = ~ 2 ( a 2,fon + a -2,£m), Q>1 m = ~ (^2, im ~ a -2, im) ■ (4-2) 

Using the af^ and it is possible to construct, in real space, A E (n, r) and A B (h, r) 

A E (n,r) = J2Ni 1 afJ l Y em (n), A B (n,r) = £ Nf 1 afiY im {n), (4.3) 

im im 

where Ni = J(£ — 2)1/ (£ + 2)!. The fluctuations A E (n, r) and A B (n, r) are the analog of 
the brightness perturbations for the intensity and for the circular polarization since they 
both admit a regular expansion in terms of ordinary (i.e. spin 0) spherical harmonics. It is 
practical to define a set of ladder operators raising the spin- weight of a given function. When 
fluctuations of different spin are treated on the sphere, a symmetry 0(4) arises naturally 
(see, for instance, [37]). Three (out of six) generators of 0(4) correspond to the generators 
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of the three-dimensional rotations, while the remaining three generators (commuting with 
the "orbital" angular momentum) are used to define the ladder operators which either raise 
or lower the spin-weight of a given function on the sphere. In the form which is suitable for 
the present calculation the mentioned operators read 



Kl(n) 



— (sin$) 



±s 



d*± 



sin 1/ 



(sin i?) 



(4.4) 



By using Eq. fj4.4j) the E-mode polarization and the B-mode polarization are given by: 



v E (n,rJ 



--{K^(n)[K {2 \n)A + (n,T)} + K^ 1 \h)[K^ 2 \h)A.(h,T)]}, (4.5) 



A B (n,r) = ^{K W (n)[^ 2) (n)A + (n,r)] - K^ 1 \h)[K^ 2 \h)A-(h,r)]}. (4.6) 

In the vanilla ACDM scenario assumed in Eqs. (12.7j) . (12. 8ft and ( 12. 9p the tensor mode power 
spectrum does not contribute; thus Ab (n,r) = implying, from Eq. (14. 5ft . that 

A E (n,r) = -^[(l- / i 2 )A P (n,r)], (4.7) 

where /x = cosi/ 1 . Recalling that A P (n, r) obeys (in Fourier space, Eq. (I2.2p ) can be 
written as 



,( E ) 



A 7 / 



dnY* m (h) 



d 3 kdl 



/i 2 ) A P (fc,/x,r) 



(4.? 



,m (2tt) 3 /2 

As anticipated, the solution of the equation for Ap(fc, /x, t) contains two contributions: one 
stemming from Sp(k, r) (i.e. the standard adiabatic term) and one proportional to Ayi- By 
solving Eq. ( 12.2ft with the line of sight integration we obtain: 

3 r T o 

Ap(fc,/x,To) = -(l-^ 2 ) / dTK,{T)S v {k,T)e-^ x 

4 Jo 



1 /" r o 
+-^ c (l-/x 2 ) / /C(r)A vl (fc,r 

2 Jo 



(4.9) 



where, as usual, x = k(r — r). To estimate Eq. ( 14. 9 p the tight coupling expansion can be 
used; the quadrupole of the intensity can then be related to the dipole of the intensity (but 
computed to zeroth order in the tight coupling expansion). As pointed out long ago this 
approximation is highly inaccurate [H]. Instead of using the dipole it is more accurate to 
obtain an equation for Sp(k,r), solve the obtained equation and then approximate the line 
of sight integral (see, e.g. [4"T| l4"2]). 

Since Sp = Aj 2 + Ap 2 + Ap it is useful to derive the evolution equations of Ai 2 , Ap 2 
and Ap by taking the appropriate moments of the corresponding brightness perturbations, 
i.e. respectively, Eqs. ( 12. ip . (12. 2p and (12. 3p . The result of this straightforward but lengthy 
manipulation is given by: 



a; 



PO 



-[A P2 + A H 



A 



POj 



-fcApi + if e e'A Y1 , 



A'i2 + e' 
A' P2 + 4 



9 1 

— A K (Ap + A P2 

10 10 V 

' 9 1 

— A P2 - — (A P0 + A I2 ) 



10 



10 



--kA 13 + hA n + y e e'A Y1 , 

3 2 
--/cA P3 + -fcA P1 + if c e'A vl . 



(4.10) 
(4.11) 
(4.12) 
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Summing up Eqs. (|4.10|) . ( 14. lip and ( 14.120 the wanted equation is given by 

S'p + ^S P = ~k(A n + A P1 ) - ~A;(A I3 + A P3 ) + 7 ~ie'f e A Y1 . (4.13) 

Note that, in the limit f c (u>) —> Eq. ( 14. 13ft reproduces the analog equation derived in 
|41j . The solution of Eq. ( 14.130 can be obtained with the line of sight method and the 
dipole of the intensity of the radiation field can be estimated to lowest order in the tight- 
coupling approximation. In the standard case this procedure gives a rather good numerical 
agreement when estimating the polarization autocorrelations as it will be clear from the 
following considerations. The E-mode is the sum of two uncorrelated contributions 

«S = C + C (4-14) 
whose associated power spectra are given by 

C (EE) = B (EE) +B CW) ! (415) 

er ~E( lasso, *r ~E( iCi 2 )- (4.i6) 

m=—l m=—l 

From the solution of Eq. (14. 13j) in the tight-coupling approximation [H] (see also [12]), the 
V-mode contribution to the full power spectrum is 

£f v) = (4tt) f 2 c \ 2 £(£ + 1)(£- 1)(£ -2) J ^p v (k) 3 -M- cos 2 (k c sv r.) e~ 2 %, (4.17) 

where A = 13/2 is a numerical factor here estimated to first order in the tight-coupling 
approximation. The diffusive damping scale £^ does depend upon the pivotal parameters of 
the ACDM scenario: 

£» = k»D A (z*) = 2240 d A (sQ ^o^-o.ii ^ (41g) 



y/r R * + 1 - y/r^ 



where, following the customary notation, = f2bo and com = Qmo (with Qyi = 
Qco + Qbo). The (comoving) angular diameter distance at z* has been rescaled, in Eq. (14.180 

as 

D A (z*) = -_— d A («.). (4.19) 

V 

Always in Eq. (14.180 the quantity r^* defines ratio of the radiation and matter energy 
densities at z*, i.e. 



Pr(z*) a eq 2 _ x f z*\ 

tr* = — f-^r = — = 4.15 x 10 u M — r . 4.20 
Pm{z*) a* V10-V 
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The numerical content of Eqs. (I4.18p - fl4.20p is fully specified in terms of whose explicit 
form can be written as 



z, = 1048[1 + (1.24 x 10- 3 )^ a738 ][l + (?rf], (4-21) 

0.0783 c^ 0238 _ 0.560 

91 ~ [1 + 39.5 a; b °- 763 ]' 92 ~ 1 + 21.1 u; b L81 ' ( } 

Equations f l4.2ip - fl4.22p represent a rather handy analytical expression for values of the 
parameters close to the ACDM best fit. This kind of approach has been also used in [32] in 
a related context and earlier analyses can be found in [331 03]. To test Eqs. fl4.2ip - fl4.22p we 
can evaluate them by using the fiducial set of parameters reported in Eqs. (12.71) and (12.81) 
and corresponding to the best fit of the WMAP 7yr data alone in the light of the vanilla 
ACDM paradigm (which is the one also assumed here). For instance Eqs. (I4.2ip and ( I4.22p 
imply z* = 1090.77 which is within the error bars of [T6"||2~T]. i.e. z* = 1090.79^0 92- Similarly 
Eq. (I4.19P implies Dx(z^) = 3.3397 /H = 14101.6 Mpc again within the error bars of the 
WMAP 7yr data implying D A (z*) = 141161^3 M P C - Equations f l4TT8l) -f T4T22l ) depend only 
upon the parameters of the ACDM scenario; Eq. (I4.17P can now be rewritten by setting 
x = w £, as 



b7 Y) = (4tt) f c X 2 £(£ + W-W -2) Ay (£) 



CO //7/1 _2-L-li) 2 

— (£w) n ^j 2 A£w) cos 2 h v £w)e ? . (4.23) 
1 w 



In the limit £ > 1 and w > 1 the spherical Bessel functions can be approximated as [33, 3 

lit + mwi) * £{£ + l ff [ ^ W 2 4) } c 1 \ (4.24) 

£ z w\Jw z — 1 2. w\Jw z — 1 



where the argument of the cosine, i.e. (3(w,£) = £\fw 2 — 1 — £arccos (l/u>) — |, leads to 
rapidly oscillating contributions averaging to 1/2 in the final expression. Similar techniques 
are employed in the semi-analytical discussion of the temperature autocorrelations [42J (see 
also [33 |36j). After some algebra B^ W ^ becomes 



l(l+l) p (W) _,2 f 2 A (o.o ,n v -ifh 



^ B? w > = X 2 fiAv(e + £v) nv -\f ) [Ve(£ D ,n v ) + Q e (£ D ,n v )], (4.25) 



n v — 1 



where 



poo w n v -7 _ 2 lL w 2 

V e (£ D} n v )= dw^==e ^ , (4.26) 

poo w n v -7 

Qe(£D,n v )= dw -—= cos (2 lv £w)e £ d . (4.27) 
Ji yjvr - 1 
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The integrals appearing in Eqs. ( I4.26P and (I4.27P shall be performed numerically by chang- 
ing, for instance, the integration variable as w = \fy 2 + 1: 

V e (£v,n v ) = - e - 2 ^/^ 2 f°° (y 2 + 1)K~8)/ 2 e -^l^?v\ ( 4 . 28 ) 
2 Jo 

Qe{£~D, n v ) = i e- 2 ^l^ 2 [°° (y 2 + l)^/ 2 cos (2£ lv Jf + l) e - 2 ^)V. (4.29) 



2 

Interestingly enough, the final result can be parametrized as 



+ D 



^£f v) = Myy(£ + ^v) n "" 1 (a v + by cos [2 lv {£ + £ v )])e- 2( ^ D)2 , (4.30) 



2tt 

AT vv = A 2 / e 2 M^(^)^ 1 T 7 2 (4.31) 
a v = 0.54, by = 0.22, £ v = 65, (4.32) 

where j v = (t^/tq)^ ~ O.OG/v^- The obtained results should be compared with the 
corresponding expression for the B\ power spectrum. Modulo the difference in the fiducial 
set of parameters (which coincide here with the ones of Eqs. (12.71) and (12. 8ft ) the result can 
be read off directly from 

£(£ + 1) 



2tt 



#f E) =Af { EE ) (£ + £E) ns+1 {aE-b E cos[2 lA (£ + £ E )]^- 2 ^ 2 , (4.33) 



A/ee = 4.453 x 10- 4 ( ^Y S ~ ( ,/* (/iK) 2 , (4.34) 



k p J V2.43 x 10 



a E = 0.194, b E = 0.110, £ E = 65. (4.35) 

The fact that £ E ~ £y is a consequence of the use of the same asymptotic expression of the 
spherical Bessel functions in both sets of integrals. In Eq. ( 14 . 3 3 1) 7a is nothing but ti/£a 
where £a is the acoustic monopole, i.e. 



UO 3 , 



In 



(4.36) 



i+v'm.-Rb, 

where Rb{z*) is the baryon-photon ratio, c s b(-2*) is baryon-photon sound speed and 



CbOO = 1 1 #b(z*) = 7- = 30.36 u h (—), (4.37) 



3[l + i? b (^)] 4 Pi 

dr 

r B (z.)=l drc sh (T)=l (4.38) 

/3[i? b (r) + 1] 



The parametrization of Eq. ( I4.36P implies that £a = 301.674 which is consistent with the 

WMAP 7yr data giving £& — 302. 5T In Fig. [2] the various approximations discussed so 

far in this section are summarized. 
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Figure 2: The semi-analytical results of Eqs. (O5]l - (j05] l and of Eqs. ( fO0j) -( fO2j) are 
illustrated, respectively, in the left and in the right panels. Since Ay is the quantity we 
ought to bound, in the plot at the right arbitrary units Ay = 1 have been employed just for 
illustration. 



In the plot at the left of Fig. [2] the full line illustrates the EE correlation computed 
numerically using the best fit parameters derived from the WMAP 7yr data (see Eqs. (12.71) 
and (12.81) ). In the same plot the dashed line denotes the semi-analytical result of Eq. (14. 33 p . 
Finally, always in the plot at the left, the dot-dashed line is the graphical illustration of the 
result of Eq. (14.331) but with &e = 0. The maximum of the dot-dashed curve is located for 



/(E) 



1 + K + 1)4 



(4.39) 



By evaluating Eq. ( P3l) in for £ = £^ the maximum of £(£+1)3^^ / (2n) can be estimated. 
In the plot at the right, always in Fig. 121 the V-mode contribution is illustrated with the same 
logic used in the plot at the left: the semi-analytical expression of Eq. (I4.26P is compared 
with the numerical results. The maximum of £{£ + l)i3( vv V(2vr) is located for 



^£\ + (n v - 1)£ 2 D - £y 



(4.40) 



To obtain the coveted analytical bound on the V-mode contribution it is thus then sufficient 
to require that: 



£(£ + l)£( vv ) 



2tt 



(V) 



< 



£(£ + 1)£( EE ) 



2tt 



max 



(4.41) 



By enforcing the condition given by Eq. (I4.4ip we obtain: 



Or 

Ay < 0.151 x — 



(jl s _|_ j^Jis+l /4e\ ( n v-n s )/2-l 



X 



A 



n 



2.43 x 10" 



ay \\ (n v - l)"-- 1 \£ D 

/ U4116Mpc/ VnG/ VGHzJ 



1091.79, 



7 v \ 2 



(4.42) 
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for n v > 1 and 

a E 



Ay < 0.151 x ^ (2£ v ) 1 -"" (n. + 1)(^ +1 )/ 2 ^ 
Arc \ / ^* + 



2.43 x lQ-V V1091.79/ V14116Mpc 



4 e \ -(n„+l)/2 

J_\ 2 /_PaW_N-- / Bu N 2 (jM 2 (4 43) 
.797 V14116Mp C y Ug7 VGHz/ ' 1 ' ; 



for < 1. 



5 Limits on the V-mode autocorrelations 

The contribution of the circular polarization to the TT and to the EE correlations has been 
separately computed in the two previous sections. Two different sets of bounds have been 
derived from complementary considerations. It is now interesting to summarize the whole 
discussion in the light possible (direct) experimental limits on the V-mode autocorrelation. 
Current bounds on circular polarization coming from direct searches are, to the best of our 
knowledge, still the ones given in [47] (see also (HI EH]) and in (50] (see also [5U [52]). The 
measurements of [50l |5T| [52] were conducted for a typical wavelength of 6 cm (corresponding 
to v — 4.9 GHz) and used the Very Large Array radio-telescope in Socorro (New Mexico). 
Conversely the limits of [471 SSI EH] used a v = 33 GHz radiometer (corresponding to a 
wavelength of 9 mm) which used a Faraday rotator to switch between orthogonal and linear 
polarization states. Besides the difference in frequency the two experiments also probed 
different angular scales. The experiment of Ref. [47] was sensitive to pretty large angular 
scales (i.e. $ ~ O(10) deg) and obtained 

y 2?r C i ; <«imK, (5.1) 

where lmK = 10 _3 K and a\ ranges from 20 to 0.2 (if we consider the sensitivity per each 
beam patch of 7 deg). The results of [501 EB E2] hold instead for a much smaller range of 
angular scales ranging from 18 arc s to 160 arc s. The authors of [50] report their results 
for the amplitude in terms of the equivalent temperature fluctuation (for which only upper 
limits existed at that time) and in terms of a putative CMB temperature of 2.75 K. In terms 
of these quantities we can establish, within the notation of Eq. ( 15 that 



where, now, c*2 is given by 

«2 



/ l(l + l) ^(W)^ „ 

^—j-—C\ ; <a 2 mK, 



(5.2) 



= 0.605, 18 arc s < ■& < 160 arc s, (5.3) 
= 0.286, 36 arc s < d < 160 arc s, (5.4) 
= 0.173, 60 arc s < d < 160 arc s. (5.5) 
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The conversion between different notations can be performed by noticing that the two-point 
function for the V-mode in real space is given by 

(Av(rn)A v (n 2 )) = £ QL±A cf v) P £ (tf), d = n x ■ n 2 (5.6) 

and by recalling that £(£+l)C^ V) /(2tt), measures, effectively, the amplitude of the two-point 
function per logarithmic interval of multipole: 

? ?£icr=/f^cr. (5.7) 

The limits listed in Eq. 05. ip and in Eqs. 05.3j) . ( 15. 4 p and ( 15. 5 p are not so recent and we are 
therefore confronted with various potential ambiguities related not only to the angular scale 
and to the frequency channel but also to the very values of the CMB temperature which 
differ from the ones determined in the framework of the WMAP 7 data release. To cope 
with this situation the following parametrization will be adopted for the direct limits on the 
V-mode power spectrum: 



^ Q (VV) = a T 70 , T 70 = 2.725 K. (5.8) 



V 2vr 

Different values of a will correspond to different observational limits either already obtained 
or potentially interesting for the present considerations. It is actually scarcely disputable that 
the observational limits on the V-mode polarization quoted in [UJ HSJ H9] and in [501 EB E2] 
can and should be improved. Using the line of sight integration, Eq. (12.31) implies 



a iV) 



— *j J dnY e * m (n) J d 3 k A v ( k, /i,r ), (5.9) 

A v {k,ii,r ) = -yz/i H drJCMe-^AviiKT). (5.10) 
4 Jo 



The limits on small angular scales on the V-mode polarization are of the same order of 
the limits obtained for larger angular scales. The V-mode autocorrelation can be neatly 
computed in the sudden decoupling limit where the coefficient a^2 can De written as 



°22 = fk^ Sm °]lwTl I d " k r° [£j£ - l(x) " (^^wWlAvi^r)^; (5.11) 

as usual, x = k{r^ — r). In the large-scale limit (i.e., in practice for £ < 40) the angular 
power spectrum of the V-mode polarization is given by: 

= 4( J+ 1)2 f T VvikMe-iix) - (£+l) 3e+1 (x)] 2 . (5.12) 
The latter expression can be explicitly computed and the result is 
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Figure 3: The bounds on the V-mode power spectrum are illustrated in terns of the amplitude 
and of the spectral index. The starred points correspond to the bounds arising from Eqs. 
(I3.45P and (13.461) . The full lines corresponds to Eq. f)5.8p with the values of a reported in 
each legend. 



cT , = ^v(|)"" I V(/,».), (5-13) 

v ^ = (STij^^ + WTW V ^ ~ WTW VM (5 - M) 

where the functions V\(£,n v ), V2(£,n v ) and Vs(£, n v ) are given by 

"P/3 Bv-\v( P 3 I n v 

r°° 1 I 2 2 I \ 2 2 

V 1 (l,n v )= dxx^-\jl_ xl2 {x) = —^ LA ^ (5.15) 

p/3 n v)~r(p 1 \ n v 
00 . o 2 . 1 V 2 2 I 1 \ C 2 _r 2 

h+3/2( 



V 2 (£,n v )=l dxx n *- 3 Jl 3/2 (x) = -^ r ^ 1 If, (5.16) 
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oo 



dx x 



3 Jt+3/2{x)Ji-l/2(x) 




-5.17) 



As previously done, it is practical to deduce a simplified expression valid in the limit £ > 1: 



Since Ay has been independently bounded from the analysis of the TT and of the EE angular 
power spectra, the l^-mode angular power spectrum is also bounded. In Fig. [3] the bounds 
stemming from the V-mode contribution to the TT power spectrum are summarized for 
different values of the magnetic field intensity. In all four plots on the vertical axis we report 
the common logarithm of \/Ay while on the horizontal axis the corresponding spectral index 
is illustrated. 

To derive the plots of Fig. |3] it has been assumed that all the parameters of the (vanilla) 
A CDM scenario are fixed to the values of Eqs. (I2.7l) - (l2.9p . The latter statement simply 
stipulates that the V-mode polarization does not affect directly the determinations of the 
ACDM parameters at least in the first approximation. The frequency channel is part of the 
observational set-up and it is therefore fixed by the nature of the experimental apparatus. 



In Fig. |3]with the full line we report the limit on \/Ay stemming from Eq. (15. 8ft in terms of 

the corresponding value of a. The various curves are obtained by using, at the left hand side 
of Eq. (I5.8P the expression of Eq. (I5.14p appropriately averaged over the multipole range. 



Always in Fig. |3]the starred points correspond to the bound on y/Ay derived in Eqs. ( 13.451) 




The results of Fig. [3] suggest that for sufficiently large frequencies and for sufficiently small 
magnetic field intensity the bounds derived from the TT correlations are not competitive 
with direct limits. This aspect can be appreciated from the two bottom plots of Fig. [3] where 
already a value a = 1 would imply a more stringent limit on \J Ay. Notice that the allowed 
region is below the full line (if the limit of Eq. (15. 8p is considered) or below the starred 
points if the limit of Eqs. (13.451) and (13.461) is enforced. In connection with Fig. [3] there are 
three possible situations: 

• the full line could always be above the starred line: this never happens in the case of 
Fig. [3] but it would simply mean that any indirect limit is more stringent than the 



• if the full line is below the starred line the indirect limit from the TT correlation is 
always compatible with the direct searches: this always happens if the magnetic field 
is sufficiently small (see, e.g. Fig. [3] bottom right plot); 




(5.18) 





and (pSED - 



direct one; 
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Figure 4: The bounds on the V-mode power spectrum as they arise from the EE correlations. 
The bounds derived in section H] are illustrated and compared with the potentially direct 
bounds parametrized, as in Fig. [31 in terms of different values of a (see Eq. (|5.8p . 



• finally the full line may crosses the starred points: this is the most realistic situation 
in the light of the present and forthcoming direct limits on circular dichroism. 

By looking at the top right and at the bottom left plots of Fig. [3] it is then apparent that, 
depending upon the frequency of the experiment, magnetic fields B n = (9(100 nG) can be 
directly excluded for v ~ GHz and with a sensitivity a ~ 10~ 6 which would imply, in terms 
of Eq. ( 15. 81) . direct upper limits on the V-mode power spectrum O(fiK) for £ < 40. 

It should be stressed that magnetic fields larger than the nG are now indirectly excluded 
by the joined analysis of the measured TT and TE power spectra [42, 53J. In this case a 
direct measurement of the V-mode polarization allows for independent constraints on the 
magnetic field intensity. At the same time the constraints derived in [4"2| [53] refer to fully 
inhomogeneous magnetic fields at large scales. Here we are constraining the uniform mag- 
netic field affecting the scattering process. The field is clearly the same but its uniformity 
is just a consequence of the smallness of the photon wavelength in comparison with the 
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inhomogeneity scale of the field. The inhomogeneity of the magnetic fields will also induce 
a Faraday rotation term which will add a coupling between the Aq and Ay. The Faraday 
term will ultimately rotate the linear polarization of the CMB [5U |55]. In this investiga- 
tion possible contributions from Faraday mixing will be ignored but their effects should be 
taken into account for more accurate estimates by following, for instance, the semi-analytic 
methodology introduced in [53] . 

In Fig. H] the same absolute bounds illustrated in Fig. [3] are now compared with the 
bounds derived in section H] from the analysis of the EE correlations. The bounds stemming 
from the EE correlations are numerically more significant, especially for large spectral indices 
(i.e. n v > 1). As in the case of Fig. [3] sufficiently small values of the magnetic field intensity 
make the indirect bounds rather loose in comparison with direct limits. There are however 
numerical differences. From the top right and bottom left plots of Fig. H] magnetic fields 
B u = O(10 nG) can be directly excluded for v ~ GHz and with a sensitivity a ~ 10~ 6 . 
The bounds stemming from the EE correlations are therefore more stringent than the ones 
derived in the case of the TT correlations. 

It is finally appropriate to mention that the frequency range assumed in the present 
discussion is in the GHz range because previous bounds, even if loose, were set over those 
frequencies. It is however tempting to speculate that, in a far future, CMB measurements 
could be possible even below the GHz. In this case direct bounds will certainly be more 
stringent but huge foregrounds might make this speculation forlorn (see, in this connection, 
[561157]). 

6 Concluding remarks 

A primordial degree of circular dichroism, uncorrelated with the standard adiabatic mode, 
can be constrained if, prior to photon decoupling, the plasma is magnetized. The obtained 
results suggest the following considerations: 

• if a V-mode power spectrum (not correlated with the adiabatic mode) is present prior 
to matter-radiation equality both the TT and the EE power spectra are affected in a 
computable manner; 

• constraints can then be inferred on the amplitude and spectral index of the V-mode 
power spectrum; 

• improved direct experimental limits on the VV correlations could be used for setting 
a limit on the magnetic field intensity. 

For experimental devices operating in the GHz range, direct limits on the circular dichroism 
imply constraints on pre-decoupling magnetic fields in the 10 nG range. Conversely, the 
current limits on large-scale magnetic fields derived from the distortions of the TT and TE 
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correlations (in the 0.1 nG range) are compatible with current bounds on the primordial 
dichroism. Improved bounds on the V-mode polarization are not only interesting in their 
own right but they might have rewarding phenomenological implications. Direct limits on the 
V-mode power spectrum in the range (9(0.01 mK) imply limits on Ay ranging from O(10~ 8 ) 
to O(l0~ 4 ) depending on the value of the spectral index and for large angular scales, i.e. 
larger than 0(1 deg). 
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A Explicit calculation of some angular integral 

In the explicit evaluation of various correlation functions it is often required to compute 
integrals involving the product of spherical harmonics and of various powers of /i. Instead of 
computing all these integrals one by one we shall just refer to the general technique employed 
here. Consider the class of integrals of the type 

dhfifi) Y e * m (h) e-** (A.l) 

where /(//) is, for practical purposes, a polynomial in /x and where, as usual, dh = d/idip. 
The integrals of the type (lA.lj) can be evaluated by noticing that /(/x) can be replaced by 
f{id x ) where d x denotes a derivation with respect to x. Alternatively it is possible to expand 
the exponential in Rayleigh series so that 

„ oo 

J dnf(fx) F/ m (n) e-*« = £H) J (2j + 1) ., ; (.r )///„.. (A.2) 



3=0 



where the term Ul m is 



i m 

2w 



Ulm = /_ x d/ijjdkp MP^Y^in, <p). (A.3) 

If f(fi) is a polynomial in /i the integrand of Eq. flA.3[) can be simplified by performing 
directly the integration over ip and by using the well known recurrence relation for the 
Legendre polynomials: 

(j + l)P,(/i) = (2j + ViiPjQi) -jPj-M. (A.4) 

Consider, for instance, the case where /(//) = 1 + /i 2 . In this case, using Eq. (IA.40 . Eq. 
(1A.3j) becomes 

U\ m = ^An{2£ +l)5 m0 ^A(£)8 tj + B(£)5j (t-i) + C(£)8j (€ _ 2) 

where 5 a b is the Kroeneker delta function and where 

(2£ + 3)(5£ 2 -!) + {£ + \f(2£ - 1) 



(A.5) 



A(£) 



(2£+ 1) 2 (2£- l)(2£+3) 



B(£) = i{e - 1] C(£) = (* + !)(* + 2) (A Q) 

y) (2£-l)(2£ + l)(2£-3Y K} (2£ + 3){2£ + 1){2£ + 5) ' y ' } 



Inserting Eq. flAT6|) into Eq. flA~2l 



dhf(jx) Y; m (n) e"** = v /47r(2£+l)5 m0 (-z) 
where 



A(£)j^) - B(£)j e ^ 2 (x) - C(£)j e+2 (x) 

(A.7) 



A(£) = (2£ + l)A(£), B{£) = (2£ - 3)B(£), C{£) = (2£ + 5)C(£). (A.8) 
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